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Abstract 

The cascade, lambda and vee type of three-level systems are shown to be described 



> 

' by three different Hamiltonians in the SU{3) basis. We investigate the Bloch space 

, structure of each configuration by solving the corresponding Bloch equation and show 

' that at resonance, the seven-dimensional Bloch sphere 5^ is broken into two distinct 

I subspaces S'^xS^ due to the existence of a pair of quadratic constants. We also give 
a possible representation of the qutrit wave function and discuss its equivalence with 

I the three-level system. 
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The atomic system interacting with two or multiple number of lasers provides an 
unique opportunity to manipulate them because of its well-defined level structure. To 
develop a systematic theory of the coherent control mechanism, it is necessary to study 
the dynamics of the two-, three- or multi-level systems which may pave the way to the 
experimental realization of the quantum computer [1-3] and associated phenomena such 
as quantum teleportation [4], quantum cryptography [5], dense coding [6,7], quantum 
cloning [8] etc. In the quantum information theory parlance where the qubit is treated 
as the primary building block, the Bloch space representation of the two-level system is 
significant for many reasons. Most importantly, a qubit state \qB > is codified as the 
superposition of the computational basis {|0 >, |1 >} and it corresponds to various points 
on the two-dimensional unit Bloch sphere [3,9]. Thus a quantum mechanical gate, 
which indeed represents an unitary operator A, corresponds to the transformation of the 
qubit state as \q'^ >= A\qB >, defined on the surface of that sphere. A natural extension 
of the qubit is the qutrit system, where the computational basis are defined as {|0 >, |1 > 
and |2 >}, respectively. Physically a qutrit possesses much complex but richer structure 
than the ordinary qubit and there exists several suggestions that it may be implemented 
either as a three-level system [10,11], or transverse spatial modes of single photons [12], or 
polarization states of the biphoton field [13,14], respectively. In recent years considerable 
progress has been made to understand various aspects of the qutrit system treating it a 
three-level system driven by bichromatic laser fields. These studies include, the quantum 
information transfer between qutrits [15], separability among joint states of qutrits [16], 
entanglement sudden death for qubit-qutrit composite system [17], cooperative mode in 
the qutrit system with complex dipole-dipole interaction [18], teleportation between two 
unknown entangled qutrits [19], cloning the qutrit [20], possible algebraic and geometric 
structure [21-24], new quantum key distribution protocols dealt with qutrits [25,26] etc [27]. 
In spite of these progress, understanding the Bloch space of all three-level configurations 
and the wave function of the qutrit system is the necessary prerequisite to define various 
logic gates and the associated circuits, which perhaps form the basis of developing a richer 
form of quantum information theory beyond qubit. 

A systematic study of the N-lcvcl system and its connection with the SU (N) group was 
first initiated by Eberly and Hioe [28-30] and later it was investigated in detail in context 
with the three-level system [31-38]. From these studies it is revealed that the quadratic 
casimir for the SU{3) group is manifested through the existence of a set of nonlinear 
constants which gives rise to the non-trivial structure of the Bloch space, i.e., iS'^XiS^XiS'*. 
Later these constants were studied by solving the pseudo-spin equation [35] and also by the 
Floquet theory technique [37,38]. However, all three-level configurations, namely, cascade, 
lambda and vee type of system, are found to exhibit same set of nonlinear constants. The 
primary hindrance of getting different constants for different configurations is the following: 
It was pointed out by Hioe and Eberly that if the position of the intermediary energy 
level E2 is changed with respect to other two levels shown in Fig.l, then all three-level 
configurations can be expressed by similar Hamiltonian [29] . Since the proposal was made, 
it has become a widely accepted convention that the cascade, lambda and vee type of three- 
level systems are described by similar Hamiltonian with the energy hierarchy conditions, 
namely, Es > E2 > Ei, E2 > E3 > Ei and Ei > E2 > E3, respectively. Similar structure 
of the Hamiltonian leads to same Bloch space for three systems and therefore forbids the 
existence of different non-linear constants for different configurations. However, it is now 
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well-known that each configuration of the three-level system is associated with a diverse 
range of coherent phenomena and therefore must be described by the distinct Hamiltonian. 
For example, the lambda configuration exhibits the phenomena such as STIRAP [39], EIT 
[40] etc, while the vee configuration is associated with quantum jump [41], quantum zeno 
effect [42], quantum beat [43] and others, respectively. Taking the interacting fields to be 
either classical or quantized bichromatic fields, the necessity of the inequivalence between 
the lambda and vee systems is also studied [44]. In a recent work, we have proposed a 
novel scheme to write the Hamiltonians of the three-level systems in the SU (3) operator 
basis where the energy levels are arranged as > E2 > Ef {a = H, A and V) shown 
in Fig.2 [45]. In particular, we have discussed the exact solution of the semiclassical and 
quantized three-level systems and have compared the corresponding Rabi oscillations. 

The primary objective of this paper is to discuss the structure of the Bloch space of 
the three-level systems, while taking distinct Hamiltonian for different configuration. It is 
shown that at vanishing detuning frequency, the Bloch space of each three-level system is 
broken into two distinct sectors <S^x<S^ and they are different for different configurations. 
Finally we have constructed the wave function of the qutrit system and discuss its various 
properties. 

It is customary to view the three-level system as two two-level systems with one of 
their levels overlaps with each other. Thus in the dipole approximation, the Hamiltonian 
of the semiclassical cascade system can be expressed as [45] 

H- = nioJiUs + UJ2T3) + di2 • Ei(Oi) + d23 • £2(02), (1) 

where, di2 = (U_|- -|- U_) and d23 = (T_|- -|- T_) be the dipole operators representing the 
transitions 1 -H- 2 and 2 o 3 and Ei(f2i) = Ei+ -|- Ei_ and ^2(^2) = E2+ + E2- be the 
bichromatic classical electric fields with frequencies Jli and ^2, respectively. Similarly we 
have 

H^ = h{0JiV3 + 0J2n) + di3 • Ei(J)i) + d23 • E2(J^2), (2) 

for the lambda system, where dis = (V+-|-V_), d23 = (T+-|-T_) be the dipole operators 
representing transitions 1 o 3 and 2 -f-^^ 3, respectively. Finally for the vee system we have 

H^ = HuiVs + u;2U3) + di3 • Ei(Oi) + di2 • £2(^2), (3) 

where di3 = (V+ + V_) and di2 = (U+ -|- U_) characterize the transitions 10 3 and 
1 -H- 2, respectively. In Eq.(l-3), ?iloi{= —Ef),fi{u:i — 0J2){= -E'2) ^'^d ^2i= Ef) be 
the energies of the three levels of the cascade system, huji{= —Ei),fiu}2{= —E2) and 
h{uj2 + ^i)(= ^3 ) be the energies of the lambda system and 7iaji(= EY), fiuJ2{= E^) and 
—h{LO2 + 0Ji){= E^) be those of the vee system shown in Fig.2, respectively. We note that 
in all cases, unlike the conventional scheme shown Fig.l [29], the energy levels maintain 
the hierarchy, E^ > E^ > Ef {a = S, A and V). In Eqs.(l-3), we have defined the SU{3) 
shift vectors, 

r± = i(Ai±zA2), c/± = i(A6±a7), v± = ^{X4±iX5), 

T3 = \3, Us = {V3\8-X3)/2, 1/3 = (V3A8 + A3)/2, (4) 

respectively, where the Gellmann matrices are given by. 
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respectively. These matrices follow the commutation and anti-commutation relations 

[Aj, Xj] = 2ifijkXk, {Xi, Xj} = ^6ij + 2dijkXk, (6) 

where dijk and fij^ = 1,2, ..,8) represent completely symmetric and completely anti- 
symmetric structure constants, respectively, which characterize SU{3) group [46]. 
In the rotating wave approximation (RWA), the Hamiltonians can be written as 

= n{ni -cui- U2)U3 + n{^2 -^i- ^2)"^^ + ^(a^?73 + A^r3)+ 

hniU-^ exp(— ■jfiif) + hK2T^ exp{—iQ2t) + h.c. , (7) 
for the cascade system, 

= n{ni -COl- UJ2)V3 + ^(^^2 - Wl - W2)T3 + HA^Vs + A^r3) + 

hniV-^ exp(— iOit) + hn2T^ exp(— ir22i) + h.c. , (8) 
for the lambda system, 

hKiV+ exp(— ir^it) + hK2U+ cxp(— i02t) + h.c. , (9) 

for the vee system, respectively. In Eq.(7-9), Ki {i = 1,2) be the coupling parameters and 
Af = {2(jji+oo2—^i) and Af = (a;i+2a;2— ^2) be the respective detuning frequencies. Thus 
we note that the Hamiltonian of any specific three-level configuration can be expressed by 
a subset of SU (3) matrices, namely, 



Cascade system : 
Lambda system : 
Vee system : 



U±,U3,T±,T3, i.e., {Ai, A2, A3, Ae, Ay, Ag}, 
T±,T3,V±,V3, i.e., {Ai, A2, A3, A4, A5, As}, 
U±,U3, V±, V3, i.e., {A3, A4, A5, Ae, A7, As}, 



(10) 



respectively [45], rather than all eight matrices [28-30]. Given with three well-defined 
models, we proceed to develop the Bloch equation of all three-level systems and discuss 
the constraints involving their solutions. 

Let the solution of the Schrodinger equation of a generic three-level system described 
by the Hamiltonian Eq.(7-9) is given by. 



^%t) = CS{t) |0) + cm |1) + cm |2) , 



(11) 
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where C(^{t), Cf (t) and C2{t) be the normahzed amplitudes with basis states |0), |1) and 
1 2), respectively. The exact evaluation of the probability amplitudes enables us to calculate 
the density matrix of any system given by 



(12) 



To start with, we consider the dressed wave function of the cascade system obtained by a 
unitary transformation 

= C/i(t)*", (13) 

where the unitary operator is given by 



(14) 



The corresponding time-independent Hamiltonian in Eq.(7) of the cascade system is given 

by 

ij2(o) = {-huliis + ulH-{t)u^) 



h{Af + 2A^) 

Hki 





?lKl 



■i^(2Af + Af) 



(15) 



Similarly, the unitary operator of the lambda system described by the Hamiltonian 
Eq.(8) is. 



UAit) = exp[--i{2n2 - ni)n + (2J^i - n2)V3)t], 



and the corresponding time-independent Hamiltonian is given by. 



H^(0) 



hK2 
flKi 



\h{A^-2A^) 











lh{A^ - 2A^) 



(13) 



(17) 



Also the Hamiltonian of the vee system in Eq.(6) can be made time-independent by 
using the transformation operator 



Uv{t) = ea;p[-|((2Q2 - ^i)Uz + (20i - J^2)^3)t], 
and corresponding Hamiltonian is, 





\n{2AX - 





A^] 



ln{2AX - 

hK2 



aX) 



-\n{AX + AX) 



(18) 



(19) 



Thus we have three distinct Hamiltonians with different non-vanishing entries for three 
different configurations. 

To obtain the Bloch equation, we define the generic SU{?)) Bloch vectors 



St{t) = Tr[p'^{t)M, 



(20) 
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where p'^ be the density matrix of the three-level systems given by Eq.(12) which satisfies 
the Lioville equation, namely, 

f = i[/,i/«(0)l. pi). 
Prom Eq.(20), the density matrix written in terms of Bloch vector is, 



(22) 



Substituting Eq.(22) and the Hamiltonians given by Eqs.(15), (17) and (19) in Eq.(21), 
we obtain the Bloch equation 

AQa 

^=M55J, (23) 

where M^^ be the eight dimensional anti-asymmetric matrix. For the Hamiltonian of the 
cascade system given by Eq.(15), the matrix Mf^ reads 
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Similarly for Eq.(17), the matrix of the lambda system is given by 
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and finally from Eq.(19) for the vee system we have, 
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respectively. 

The algebraic structure of the SU (3) group allows the existence of a set of quadratic 
casimirs which will appear in form of the quadratic constants. However, the searching of 
the exact tuple of the Bloch vectors forming such constants is quite cumbersome because 
of the large number of such combinations. More specifically, for eight Bloch vectors, we 
have to search j^~^{= 8, 28, 56, 70, 1) number of combinations forming such tuples 
for n = 1, 2, 3, .., 8, respectively. We have developed a Mathematica program to search the 
exact combination of the Bloch vectors obtained by solving the Bloch equation in Eq.(23) 
and after a rigorous search obtain the following results for n = 3 and ra = 5 only: 

For the cascade system at resonance (Af = = Af), the Bloch space is split into two 
parts; one being of two sphere S^, 

Sf{t) + Sf{t) + Sf{t) = Sf{0) + Sf{0) + Sf{0), (27a) 
and other is the four sphere »S^, 

Sf{t) + Sf{t) + Sf{t) + Sf{t) + sfit) = 

sfiO) + sfiO) + Sf{0) + Sf{0) + Sf{0), (276) 

respectively, where 5*^(0) be the constants at t = which are to be evaluated in terms of 
the probability amplitudes. By noting the fact that the density matrix can be written as 
p^{t) = U\t)p-{Q)U{t), Eq.(20) becomes, 

5f(0) = rr[p=(0)Ai]. (28) 

Plucking back Eqs.(5) and (12) into Eq.(28), different values of the constants 5^(0) can 
be evaluated in terms of probability amplitudes Cf^2,3(0) at t = and Eqs.(27) become 

Sf + Sf + Sf = 4|Cf (0)|2|Cf (0)|2 + 4|Cf (0)|2|Cf (0)|2 (29a) 

and 

sf + sf + sf + sf + sf = ^(|cf (o)|2 + |c|(o)|2 + \cm\^f 

-3|Ci=(0)HC2=(0)p - 3|C2=(0)HC3=(0)p, (296) 

respectively. Similarly at resonance (Aj^ = = Ag ), for the lambda system we have two 
Bloch spheres, 

sf + sf + sf = 4|(7f (0)p|C7f (0)|2 + 4\Ci{0)f\Cf{0)f (30a) 

and 

sf + ^3^' + ^5^' + ^6^' + ^8^' = lilCHof + 1^2^0)1^ + \CH0f? 

-3|C7i^(0)|2|C2^(0)|2 - 3|C73^(0)|2|Ci^(0)|2. (306) 
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Also for the vee system at resonance (Aj^ = = Ag'), we have 

Sf + + Sf = 4|C7r(0)p|C3^(0)p + 4|C73^(0)p|Cr(0)p (31a) 



and 



sx' + ^3^' + sf + sf + sf = ^(icr(o)p + \cYm' + icYmr 



-3|C2" (0)nC3^ (O)l^ - 3|C3^ (0)hC/ (0)|^ (315) 

respectively. Thus we note that the 3- and 5-tuple of the Bloch vectors form two quadratic 

constants which arc different for different three-level systems. From Eqs. (29-31) it is 
further evident that at resonance, the Bloch space 5'^ is broken into two-subspaces xS^, 
which indicates that each configuration has distinct norm. It is instructive to note that at 
off resonance (Af / / A2), the solutions of the Bloch equation of all three configurations 
satisfy, 

sf + sf + sf + 5f + sf + sf + sf + sf = ^(|Cf(0)|2 + 1^2^(0) |2 + |C^(0)|2)2. (32) 

This shows that the normalization condition |Cf (O)p -h |Cf (O)p -|- |Cf (O)p = 1 gives that 
the radius-squared of the seven-dimensional Bloch sphere S"^ will be | regardless of the 
configurations. 

Finally we consider the wave function of the qutrit system. A convenient parametriza- 
tion of the normalized qutrit wave function is 

\qT >= cos -^-10 > -l-sin — sin — sin —e^'^\l > + 

(sm — cos — +ism — sm — cos — )|2 >, (33) 

where the three coordinates 6q, 9i and 62 lie between to vr with arbitrary value of the 
phase angle (p. From Eq.(33) it is evident that, the pure qutrit states |0 >, |1 > and 
|2 > can be obtained by choosing the coordinates (^O) ^i; ^2) to be (0, ^1, ^2)) (ti"; t^, t^) and 
(7r,0, ^2)) respectively. This indicates that for a qutrit system, the state |0 > may be any 
point on the (^1, 62) plane, |1 > corresponds to a well-defined point and |3 > lies on a line 
for any value of 63, respectively. In other words, in the three-level scenario, the degeneracy 
of the |0 > and |2 > states shows that the Bloch space structure of the qutrit is highly 
non-trivial. This is in contrast with the two-level system representing qubit, where the 
basis states |0 > and |1 > correspond to the north and south poles of the Bloch sphere, 
respectively [3,9]. Using Eq.(12), various elements of the density matrix corresponding to 
the qutrit wave function in Eq.(33) are given by 

11 2 ^1 
Pt =cos — 

22 . 2 ^0 . 2 ^1 ■ 2 ^2 
Pt = sm — sm — sm — 
2 2 2 
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1 ^0 

= - (3 + cos 9i + cos 02 — cos 9i cos 02 ) sin^ — 

12 21 * 1 7fli ■ /I . ^1 . ^2 

p^^ = = -e*"? sm sm y sm — 

23 32* ■ 2 ^0 . ^1 . . 9i O2. . 02 

Pt = Pt = e ^ sm — sm — I cos \- 1 sm — cos — ) sm — 

'^^ 2 2^ 2 2 2 ^2 

13 31* 1 • /) / ^1 , ■ • ^1 

Pt = Pt =2 ^^"^ Oq(c.os — + i sin y cos y ) (34) 

It is worth noting that this qutrit wave function is normahzed, i.e., Tr[pT\ = 1 and the 
density matrix satisfies the pure state condition, namely, p\ = pr- Finally plucking back 
pT in Eq.(20) we obtain the Bloch sphere on S"^, 

Sf + Sf + Sf + Sf + Sf + Sf + Sf + Sf = \. (35) 



This is precisely same as Eq.(32) obtained from different configurations of three-level 
system indicating the equivalence between the qutrit system with the three-level system. 

In this paper we have discussed different semiclassical three-level configurations and 
obtain the solution of corresponding Bloch equations. It is shown that, if the energies of the 
three levels follow a definite hierarchy prescription, it is possible to write the Hamiltonians 
in the SU{3) basis which are different for different configurations. At zero detuning, the 
Bloch sphere 5^ of each configuration is broken into two disjoint sectors S^xS^ due to 
the existence of different constraints which are not alike. A significant outcome of our 
approach is that we have given a possible representation of the qutrit wave function which 
can be identically represented by a three-level system. Although our formulation gives a 
possible representation of the qutrit wave function, but it is quite difficult to ascertain the 
complicated structure of the Bloch space at this stage. Apart from that, the mechanism 
of minimizing the information loss while teleporting between two or multiple number of 
qutrits, developing appropriate search algorithm, the qutrit based protocol transfer etc 
are still open issues which require further exploration of the qutrit wave function proposed 
here. Although the quantum information processing has made significant progress with 
the conventional qubit as the primitive identity, but the potential application of the qutrit 
can not be properly addressed unless we do not understand the wave function of the 
quatrit and its relation with the Bloch space structure of the three-level system in greater 
detail. More exploration along this direction may provide some interesting results which 
are within reach of the future high-Q cavity as well as nanomaterial based experiments. 
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Figure 1: Cascade, lambda and vee type three- level configurations according to the 
Eberly-Hioe scheme [29] with the energy conditions, > E2 > Ei, E2 > E3 > Ei and 
El > E3 > E2, respectively. Here the position of the middle level (level-2) with energy 
E2 is changed to generate all models. 
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Figure 2: Cascade, lambda and vee type three-level configurations according to our 
scheme [45] > > Ef (a = S, A and V). 
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